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1.0 Introduction 



The purpose of this investigation was to tost the relative 
importance of various aspects of correcting predicted values on a 
grid by incorporating information from observed vraluos at scat- 
tered data points. Grid and observation configurations v>^ere 
patterned after those routinely available over North America. 
Although investigations were limited to the univariate objective 
analysis methods, I believe the results are indicative of those 
that would be achieved in the more general case. 

Previous investigations on the error contribution of various 
steps in the objective analysis process are limited. Koehler 
(1979) separately studied the errors of a number of grid-to- 
observation and observation-to-grid interpolation (approximation) 
routines. He noted that although little attention is typically 
paid to the gr id-to-observat ion interpolation process significant 
errors may be caused by this phase of objective analysis. While 
this may be a surprise since these errors are usually small 
compared to the first-guess errors at the grid points, my results 
further demonstrated that the contribution to overall error made 
by the g r i d- to-obse r va t i on interpolation process should not bo 
ignored. This investigation complements recent v/ork by Seaman 
(193'1) regarding tlie accuracy of statistical and successive cor- 
rection schemes. His work provides expected mean squared error 
estimates for these schemes. His v/ork is very thorougli in that 
it provides estimates of the analysis error as the parameters of 
the first-guess error are varied while holding the assumed values 
constant, and vice-versa. 



1 



In Ssction 2 I derive a gonf-' rci 1 i zed expression for tli"^ 
overall error in objective analysis which leads to several 
observations. In Section 3 I describe the simulation method and 
the various options which can be easily handled. In Section 4 I 
present the results of the simulations and discuss their implica- 
tions with regard to the observations made in Section 2. 

2.0 The Form of the Error Term in Objective Analysis 

My setting for study of the objective analysis process 
assumes the following: 

(i) The true field (function) to be analyzed is H. 

(ii) H is known imperfectly at grid points through a "first- 
guess" which is in error by an amount to be denoted by g. The 
error is a normally distributed stationary random function which 
has a certain spatial correlation and standard deviation. 

(iii) H is imperfectly measured at observation points yield- 
ing values with errors o. These errors are independent and 
normally distributed with certain standard deviation. 

The nature of the errors makes it only possible to evaluate 
g at grid points, and o at observation points, although it is 
sometimes convenient to think of them as functions rather than as 
sets of errors. The objective analysis process consists of int- 
erpolation of the first-guess values from the grid to the obser- 
vation points (by a linear operator designated M) followed by 
interpolation of the difference between the observed and first- 
guess values back to the grid point (by a linear operator desig- 
nated L) as a correction to the first-guess values. Denote tlie 
error in the entire process by E, then the final approximation is 
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H + E = Hfg + L(Il + o- M(H + cj)) 

Let m(H) represent the error in the approximation of ii by M(I1), 
then M (H) = H - m(H). Rearranging and simplifying the above, 

leads to 

E = g + L(H + o- M (H) - M (g) ) 

= g + L(H +o - H + m(H) - M(g)) 

= g + L(o f m(H)) - LM(g) , 
and finally, 

E = L(o) + Lm(H) + (g - LM(g)) . (1) 

Thus the error is made up of three parts. The term L(o) i .s 
dependent on the 'function' o, which describes instrumentation 
error and is typically not controlable. It is obviously advan- 
tageous to have o small. Since the values of o are assumed 
independent and random it is desirable for L to be a smoothing 
operator. The second part, Lm(H) is within our control and the 
gr id-to-observat ion point interpolations error should be made 
small. If it is, then interpolation of the error back to the 
grid points by L is also small, assuming this smoothing operator 
is typical and does not magnify the error. The third part (g - 
LM(g)) is the error in interpolation of the first -guess error at 
the grid points to the observation locations by M, then back to 
the grid points by L. While it is possible that a certain sym- 
biosis between parts could occur, the goal is certainly for each 
interpolation process to have small errors. Ideally the operator 
L should be a loft inverse of tlie operator id, although this is 
almost certainly impossible. 

Partitioning the orroi in this way shows, fot example, tliat 
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using a better interpolation process from the grid to the obser- 
vation points should decrease the overall analysis error. In 
certain realizations, of course, the errors may tend to cancel. 
Since the three terms represent uncorrelated errors, the total 
error variance over many realizations will tend to be the sum of 
the individual variances. Thus, decreasing any one will lead to 
statistically smaller error variances. 

3.0 The Computer Simulation Methods 

In order to simulate the behavior of the overall error under 
various interpolation processes and first-guess error assump- 
tions, a modular computer program was written to give several 
options for the different processes. This made it possible to 
test a large number of combinations of methods and assumptions. 

In general terms, the process simulated consists of the 
following steps: 

(i) An underlying ma thama t ica 1 ly defined function 
describing the field to be analyzed is evaluated on a 
grid of points. 

(ii) "1' i r s t-guess" error is generated from normal random 
deviates with a pr e-spec i f i ed standard deviation and spatial 
correlation. 

(iii) "Observed values" are generated by evaluating the 
field to be analyzed at the observation points, and adding 
normally distributed uncorrelated random 'leviates to these 
values . 

(iv) The first-guess values at the observation points are 
obtained by one of several interpolation schemes. 
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(v) Based on the difference betv/een first-guess and obser- 
ved values at the observation locations, "corrected" values at 
tbie grid points are obtained. I will refer to the corrected 
values as the analysis values. 

Most of the simulations were done with two different grids 
and observation point sets. One was based on a 2.5° grid cover- 
ing 112.5° W to 32.5° W and 3 0° N to 5G° N, v;ith 117 = 13x9 grid 
points and 3^1 observation points within the grid, as sho’wn in 
Figure 1. The other was based on a 5° grid covering 125° W to 
75° W, and 25° N to 50° N, with 8 3 = 11x8 grid points, and 57 
observation points within the grid. This grid and the observ)- 
tion locations are shown in Figure 2. All the simulations used 
were univariate analysis methods on a, t wo-d i mens i on a 1 field. 

This simplification was necessary for two reasons. The first 
reason is that the generation of error with a specified spatial 
correlation required factorization of the correlation matrix into 
the product of a lower triangular matrix and its transpose. The 
correlation matrix is of order equal to the number of grid points, 
and it is not particularly v^oll conditioned. I nco r j^o ra t i on of 
multiple levels, a large grid, or correlated multiple variables 
was therefore not possible. The other reason is that statistical 
results required that numerous realizations bo simulated, thereby 
limiting the time available to do the computations. 

The underlying mathematically defined field can be any spec- 
ified function. The height field test function used is the one 
given by Koehler (1979) and also described in Wahba and 
Wendelborger (1930). The input pdrametors, 0^ (the location of 
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the longitudinal wave), A9, (amount part of the field is skov/od 
logitudinal ly) , and p (the pressure for the height field) are 
easily varied. The experiments simulated the 50D mb height 
field, using fixed or randomly varying 9^ and A9. A typical 
field of height contours generated by this function is shown in 
Figure 3. First-guess errors had a nominal standard deviation, 
rg, of 30 m. The spatial correlation function was modeled using 
exp( (-d/Cjj) ) ^, where d is distance (on the degree grid), and 
is a correlation distance, specified as 10°. I have used degree 
measure for distance rather than true distance, to maintain a 
rectangular grid of first-guess points. This resulted in a 
distortion of the distance varying with location. The observa- 
tion errors had a nominal standard deviation, r^, of Id n. The 
observation locations approx i mately correspond to the North Ame- 
rican radiosonde network within the grids being used. They are 
shpwn, along with the grids, in Figures 1 and 2. 

Tiie output consisted of mean, root-mean-square, and maximum 
errors over each data set (first-guess at grid points, first- 
guess at observation locations, observation at observation loca- 
tions, and analysis values at grid points) for each realization. 
The first and third of these mainly served as a check on the 
psuedo-random number generator (IMSL subroutines GGNSM and 
GGNNL). The output also gave summaries of the same errors over 
all realizations as well as the mean and standard deviation of 
the root-inean-square errors over the realizations. Interpolation 
processes are sometimes ill-behaved around boundaries. Since in 
the global problem this can be avoided, the effects were mini- 
mized here by tabulating error only over the interior grid 



points. Thus tha results are over 77 grid points on tlie 
grid and 54 grid points on the 5° grid. The options simulated 
for each step are described below. 

a. Gr id-to-observation point interpolation 

First-guess values at the observation points are obtained by 
interpolation from the first-guess grid values. I compared four 
schemes. Others could be easily included, however my resu''ts 
indicate it will probably not be fruitful to do so. The methods 
I have used are: 

(i) Piecewise bilinear interpolation. As v/ith any (piece- 
wise defined method, one must first determine the rectangle in 
which the evaluation point lies. Then, the evaluation is most 

O 

easily seen as translation to the square [0,1] , followed by 3 

one dimensional interpolations. This requires 8 operations, 
where an operation is defined as a multiplication or division 
followed by an addition or subtraction. Practically, the evalua- 
tion can be accomplished in 5 operations (and a couple of extra 
add i t i ons/subt r ac t ions) . In my cost analysis I have used 8 
operations; this cost is very low compared to that of other 
necessary calculations. 

(ii) Bicubic spline interpolation. I used the IM8L subrou- 

tines IBCCCU and IBCEVL. Preprocessing for the spline coeffi- 
cients on a grid requires 1 2N^xM^ + 27:M^ + 51 Ng- 1 h operitions. 

Evaluation requires 2 operations to translate to [0,1] and 5 
cubic interpolations at 9 operations each for a total of 47 
operations. The preprocessing operations involve solution of 
tridiagonal systems of equations -wliich aic amenable to vectoriza- 
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tion for pipeline computers. 

(iii) Piecewise bicubic interpolation. My implementation 
of this scheme used 2 operations for a translation to ro,31 
followed by 5 cubic interpolations, each costing 5 operations. 

In addition, a difference table was formed at a cost of several 
subtractions . 

(iv) Bessel bicubic interpolation. My implementation of 
this scheme used 2 operations for a translation to [0,3] fol- 
lowed by 5 cubic interpolations, each costing 5 operations. 
Because of default to parabolic interpolation in boundary 
regions, there were some additional tests. There were also a few 
subtractions to form the difference table. 

b. Observa tion-to-g r id point interpolation 

As in operational weather forecasting programs, the differ- 
ences between first-guess and observed values at the observation 
points are used to correct the first-guess values on ttie grid to 
obtain analysis values on the grid. I have tested twelve schemes 
for performing this correction. I will give a brief description 
of each method and refer the reader elsewhere for complete de- 
tails. The first-guess error at the observation location, P|^ = 
(01^,ci>l^), is denoted by AH|^, I< = 1 ,...,Nq. The number of grid 
points is ^0^^. I want to evaluate the approx i ma t ion at grid 
points, but will write; it in terms of a generic point, P = (9,<i)). 
Recall that the standard deviation of the first-guess errors is 
Cg, and the spatial covariance function is denoted by C(P,Q). 

An operation count has been made for each of the methods. I 
discuss briefly how various phases of the process contribute, and 
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summarize ^.he results in Table I, along v/ith soma r?[)res?nta t j /. 
numbers that arise from my simulations. 1 have described soi.ts 
schemes as local, implying that othisrs are global. In the con- 
text of global objective analysis, all the schemes 1 consider are 
local; the schemes which are global for my simulation are ]ess 
local than the ones I refer to as local. 

(i) Optimum interpolation (01). This scheme was introdu"'’'d 
to the meteorological literature by Gandin (19o3) and has re- 
ceived widespread attention in recent years, e.g. see Bergman 
(1979) and Lorenc (1981). The method in its proper form requires 
that the spatial covariance function of the first-guess errors 
and the standard deviation of the observation error be known. 
Since these are known for this simulation, I have used their 
properties. I have implemented the scheme as described in Franke 
and Gordon (1983), viewing the approximation as a linear combina- 
tion of the covariance functions associated v/ith the observation 
points. Thus we have 
% 

AH(P) = E a,^C(P,P,^) , 
k = l 

where C(P,Q) is as noted above. The a|^ are determined from the 
system of equations 




where AH^ is the difference between the first-guess and observed 
values at the i*"*^ observation point, r^ is the standard deviation 
of the observation error, and 6^j is the Kronecker delta. 

The cost of (01) consists of a preprocessing phase that 
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includes the generation and solution of the systoiri of equations, 
followed by evaluation of the analysis at the grid points. For 
Nq observations, preprocessing is at a cost of Nq([s'q + 1)/2 func- 
tion evaluations to generate the coefficient matrix and 
Nq)/ 6 operations plus square roots to perform Cholesky decom- 
position and solution of the system of equations for the a|^. 
Evaluation costs covariance function evaluations and opera 
tions to form the linear combination representing the value of 
the correction at each grid point, 

(ii) Local optimum interpolation. In my version of this 

scheme, nominally only points within the surrounding 10° square 
are used; if fewer than 4 observations are available, the squar 
is expanded to 15° and so on, by 2.5° increments in each direc- 
tion until at least 4 observations are availabl.-, Ttie costs of 
the search were not assessed. For each grid value correction, a 
system of equations must be formed and solved, and the correspon 
ding correction computed. With n observations being used the 
expressions given for 01 above apply with n replacing This 

process was performed for each grid point, making the total cost 
the sum of these costs over all grid points. 

(iii) Global Barnes' method. This type of scheme is des- 
cribed by Barnes (1973) and others. My scheme used the known 
correlation functions as the weights for the first pass. Thus, 
the appr ox i ma t ion is 

o o 

AH.^l(P) = Z W^(P)AH,^/ Z W|^(P) , 
k=l k=l 

where W|^ = exp ( (- ( I I P-Pj^ I 1 /c^) ) , and AH|^ is as before. For the 
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second pass the correction has t!ie same £orm, but AH |., is r<? pi a cod 
by AH|^ , the difference between the corrected first-guess and 
the observations. The quantity c^ is replaced by for the 

second pass. The total correction at the grid points is then the 
sum of the two corrections. For each grid point the cost of this 
method is weight function evaluations per pass and ^JQ^■1 opera- 
tions per pass. In addition a separate interpolation from the 
grid points to the observation points is required before the 
second pass. This type of scheme has been defined and studied in 
a different context, without a change of weight functions between 
iterations, by Foley and Nielson (1930). 

(iv) Local Barnes' method. The same localization process 

as used for the local 01 scheme (ii) was used here. As for the 
global version, two passes were used. Hence the cost for an 
evaluation at a grid point with n neighboring points is the same 
expression as in the global scheme, but with n replacing In 

addition, there was the search cost to determine the nearby 
observations, which was not assessed. Costs of an interpolation 
from the grid points to the observation points between passes was 
included. 

(v) Statistical interpolation (c^ = 14°). In practical 
applications of 01 the error correlations and standard deviations 
cannot be modeled precisely. This iias lead to the use of the 
name "statistical interpolation". Computationally the method is 
identical to the 01 scheme (i). Here the only difference is the 
substitution of an inexact correlation distance, c^^ = 14. The 
algorithm and costs are identical. 

(vi) Statistical interpolation (c^j = 7°). Again this is 
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identical to (i) except that tlie inexact \/alue substituted for c,, 
is 7 . 

(vii) Statistical interpolation (damped cosine correlation 
function). Once more this scheme is computationally identical to 
(i) except that the correlation function used is of the form 
exp( (- ( I I P-Q I I /c^) cos ( ( I I P-Q I I /c^) ( /2)). I used the value c^, 

= 10 . 

(viii) Thin plate splines. This method is described by 
Wahba and Wendlelberger (1980) and others. The approximating 
function used by the scheme is 

% 

H(P) = Z A|^B(P,P,^) + a0 + bcii + c , 
k = l 

where the basis function B(P,Q) = | |P-Q| l^logl I P-Q 1 I . The A]^ 
and a, b, and c, are obtained by solving the system of equations 

2 

E A,^(B(Pi,Pj)+XN,^rQ ^ j ) +a0 ^ + b<i ^ +c = A H ^ ,i = l,...,dQ 

j = l 

^o 

Z Aj9j = 0 

j = l 

^o 

E Aj4)j = 3 

j = l 

^o 

E Aj = 3 . 

j = l 

In the above, A is a smoothing parameter. The smoothing parame- 
ter was chosen on tiie basis of a few trials with no attempt to 
optimize its choice for a particular data set, as can be done. 
Wendelberger (1931) describes a program that will automatically 
choose A (and m as well, see next method), but I liave not tested 
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it yet. This system of equations is symmetric, but no^ positi-u 
definite. I have used standard L-U decomposition routines to 
solve the system. Methods for symmetric indefinite systems 
about half as many operations, however I observed greater numeri- 
cal stability using the general decomposition process. Tfiere arc 
Mo(Mo+ 1)/2 basis function evaluations, and solution of the system 
of equations requires (Mq+ 3) (N^ + iSN^+S ) /3 + (Nq+ 3)^ operations. 

Unlike symmetric positive definite systems, solution of these 
equations requires searching for a pivot and pivoting. Evalua- 
tion at each grid point then requires basis function evalua- 
tions and Nq+ 2 operations to form the sum. 

(ix) Laplacian smoothing spline (m=3). This scheme is also 

described by Wahba and Wondelberg er ( 198 ?,), and is one of those 
available in the program by Wendelberger (1981). The thin plate 
spline method is a member of this family (with m = 2), but also has 
the "thin plate" interpretation. The reason for inclusion of 
this method is that the results of V»aliba and Wendelberger indi- 
cate that pressure height fields are better approximated using 
values of m = 3 or 4. I will not describe the method fully. It 
requires evaluation of Mq(iMq+ 1)/2 basis functions and multi- 

plications to set up the system of Nq+ 6 equations to be solved. 
Then Mq+ 5 operations would be required for evaluation at each 
grid point, along with tiu; evaluation of basis functions. 

(x) Franke/Gordon. This scheme was suggested by Franke and 
Gordon (1933) as one which is an explicit scheme, similar to 
Barnes* method, but which when iterated converges to the 01 
interpolant. Three iterations, with the parameter = .85||Ml| 
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(in the notation of that report) were performed. The cost in 
operations is 2 Nq(Nq+ 1) plus for each grid point. The number 

of weight function evaluations is 2 Nq plus 3 Nq for each grid 
point. 

(xi) Pseudo-Barnes' metliod. Ttiis method was described in 
Franke and Gordon (1983) and was at that time mistaken for Barnes' 
method. It differs in that the error at the second iteration is 
Barnes' approximation evaluated at the observation point minus 

the first-guess error, rather than the the corrected first-guess 
at the grid point interpolated to the observation point minus the 
observed value. The cost of this algorithm is evaluation of 
weight functions plus 2 Nq for each grid point. It requires 
No(No+l) operations, plus 2(N^+1) for each grid point. 

(xii) Local pseudo-Barnes' method. This is a local version 
of (xi), using the same ''nearby" observation points as (ii) and 
(iv). A grid point with n nearby observation points requires 
evaluation of n +2n basis functions and n +3n+2 operations. 

4.0. Results 

The simulation program described in the previous section was 
run for a substantial number of different options. Each run 
consisted of 100 realizations of a test field each containing 
associated first-guess and observation errors. Table 2 gives the 
assumed parameter values for the various cases. Not all combina- 
tions of gr id-to-observation point and observation-to-gr id point 
interpolation schemes were used in every case. The tables detail 
the complete results and the entries indicate which combinations 
were computed. Each combination in a given table (3-14) 
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corresponds Lo the same set of realizations, but different tables 
depend on different realizations. 

This investigation was designed to determine the influence 
of the g r i d- to-obse rva t ion point interpolation scheme. This 
influence is seen by noting changes in error for particular 
observe t i on- to-g r id point interpolation scheme as the grid-to- 
observation point interpolation scheme is varied. The rows of 
Tables 3-14 give this information. The bicubic spline 
interpolation produced significant improvement over pieccv/ise 
bilinear interpolation. This verifies the smaller magnitude of 
the term Lm(H) in the error expression given by (1) for the 
spline method. For 2.5° grids the errors were no smaller for 
spline interpolation than for piecewise bicubic or Bessel bicubic 
interpolation. Evidently the grid spacing was small enougli (for 
the test function used) that the interpolation error was not 
significant. Spline interpolation did show an improvement over 
piecewise bicubic and Bessel bicubic interpolation on the 5° 
grid. Spline interpolation and the cubic interpolation methods 
showed even greater improvement over piecewise linear interpola- 
tion on the 5° grid than on the 2.5° grid. Interestingly the 
first-guess errors at the observation points had greater rms 
values for cubic interpolation than they did for linear interpo- 
lation. This occurs because linear interpolation inhvercntly has 
greater smoothing. 

'^ost of the useful information given in Tables 3-14 can be 
more easily obtained from plots of the salient values. Figures 
4-8 give plots of skill vs. cost of the algorithm in thousands of 
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operations per analysis. Here "skill" is defined to be 1 - 
rmsa/tQ, where rmsa is the rms error in the analysis values. The 
skill with respect to bilinear and bicubic ititerpolation aro each 
indicated, connected with a straight line to delineate the extent 
between the two. The results for only one of the statistical 
schemes, (vi), has been plotted since the others were nearly 
identical. For these purposes I counted an evaluation of a basis, 
weight, square root, or covariance function as 10 operations. 

The plots reveal that the statistical schemes, local vOl, and thin 
plate splines all had close to the same accuracy and all were 
slightly less accurate than 01. The Barnes' schemes, the 
Franke/Gordon scheme, and Laplacian smoothing splines were least 
accurate. The poor performance of the Laplacian smoothing 
splines here, in contrast to the better performance obtained by 
Wahba and Wendelberger (1930) is probably due to the scheme being 
applied to the first-guess error function rather than to the 
underlying true height field. The degradation in tiie performance 
of the less than optimal statistical schemes is perhaps less 
drastic than one might expect. It does appear tiiat it v/as better 
to underestimate the correlation distance than to overestimate 
it. 

Figure 9 shows plots of the rms errors in the analysis 
values as a function of first-guess errors. The improvement in 
the Barnes' scheme as the first-guess errors decrease was rapid. 
The scheme gave results nearly as good as 01, the statistical 
schemes, and thin plate splines. This occurred because the 
principal problem became smoothing observation errors as the 
first-guess errors tended to zero. Figure IK shows plots of the 



rms errors in the analysis values as a function of observation 
errors. As observation errors go to zero the importance of 
modelling the first-guess error was more important than 
smoothing. Thus 01, the statistical schemes, and thin plate 
splines improved tiie most, while both Barnes' sctiemcs improved 
little. Figures 11-13 show the rms errors in the analysis values 
when incorrect variances were specified for the interpolation 
routines. Methods not using these values were naturally unaf- 
fected so that changes in tiie rms errors in the analysis values 
for these methods only reflect the variability of the (different) 
realizations used in the various cases. The plots show that the 
use of incorrect values for the first-guess and observation error- 
variances did not drastically affect the accurr4cy of the statist- 
ical methods. The interested reader is referred to Seaman (1983) 
for more extensive tests of the effects of incorrect [rarametor 
specification on the performance of statistical interpolation 
methods . 

One of the attractive features of the statistical scliemiis is 
that they afford a calculation for the estimated mean squared 
error. These estimates do not depend on any (><irlicular rcvalizi- 
tion, so they were not incorporated into the process. However, I 
did compute them as a side calculation for my grids an'! observa- 
tion points. The results of these calculations are tabulated for 
the 2.5° grid, along witti the empirical rms errors obtained 
during the simulations. Table 15 shows that the estimates given 
by 01 w’ere quite good; the estimated and empirical errors varied 
only a few percent. They also were accurate for local 01, as 
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they should be. On the other hand, the slight degradation in 
performance of statistical methods when incorrect correlations or 
variances were specified did not carry over to the error esti- 
mates. In fact the schemes that have their performance degraded 
the most (in this case, using too long a correlation distance) 
showed a decrease in the estimated error variance. Conversely, 
shortening the correlation distance in the statistical method 
increased the error estimate as well as the empirical error 
obtained, although the empirical error is underestimated. This 
indicates that one must not put too much faith in the error 
estimates when the actual covariance structure is not known, as 
in practice. It appears one could obtain just about any error 
estimate v/ished simply by specifying unrealistic parameters Cor 
the covariance structure. 

The principal results of this study were as follows. The 
decomposition of the error into independent components in (1) 
identified possible ways to decrease the analysis error. This 
lead to the results showing the contribution of the grid-to- 
observation interpolation process, the necessity of smoothing in 
the observation-to-gr id interpolation process, along with accuracy. 
The simulations provided confirmation of the above and yielded 
information concerning the sensitivity of statistical interpola- 
tion schemes to inexact parameter specification. The error esti- 
mates provided by statistical schemes were shown to be sensitive 
to inexact parameter specification. 
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seudo- 2 i'1qN^ (Mq + 1) (Mq+ 1) 2i'JQ(NlQNi — 9990, 9720 

B.-»rru's’ ir^S2'1, 19281 



Tablelf 


■^g 


''o 


®o 


0 


grid 


no lie 3^ 




3 


30 


10 


100 


0 


13x9,2.5° 






4 


30 


10 


100 


13.775 


13x9,2.5° 






5 


20 


10 


100 


0 


13x9,2.5° 






5 


30 


5 


100 


0 


13x0,2.5° 






7 


30 


10 


random^ 


random^ 


13x9,2.5° 






8 


30 


10 


random'^ 


random^ 


13x9,2.5° 


o 

II 


30 


9 


20 


10 


random^ 


random*^ 


13x9,2.5° 


II 


20 


10 


30 


10 


random*^ 


random^ 


13x9,2.5° 


(-^o^i = 


5 


11 


30 


5 


random^ 


random^ 


13x9,2.5° 


(>^o)i = 


10 


12 


5 


10 


random*^ 


random^ 


13x9,2.5° 






13 


30 


0 


random^ 


random*^ 


13x9,2.5° 






14 


30 


10 


random^ 


random^ 


11x8,5° 












Table 2 








^ ^o 

b A9 


un i f ormly 
uni f ormly 


distributed in 
distributed in 


(-82. 5°, 112. 5°) 
(-15°, 15°) 







c Tne statistical interpolation routines were given 
incorrect variances, as indicated 
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F = 50C, lheta= 100, Celth = 0 Entries: RMSE 

rg = 3C, ro = 10 Mean 

Number cf realizations = 100 

13X9 grid of 2.5 degrees, 36 observation points 



Grid-to-cbs : 


FW linear 


Bicub Spl 


PW Bicub 


Obs-to-grid 








Opt Interp 
(Cd = 10) 


6.64 

6.53 (1. 18) 


6 .09 

5 .97 (1.20) 


6.09 

5.98 (1. 19) 


Local 01 
(Cd= 10) 


7.09 

6.99 (1. 19) 


6.53 

6.42(1.19) 


6.54 

6.44 (1. 1 9) 


Barnes' 
2- Pass 


9.27 

9.08 (1.03) 


8.87 

8.68 (1.82) 


8.87 

8.68 (1.82) 


Barnes ' 
(Lccal) 


8.4 2 

8.27 ( 1. 57) 


7 ,95 

7.79 (1.56) 


7.96 

7.80 (1.56) 


Stat Interp 
(Cd = 14) 


7.28 

7.23( 1. 22) 


6 .78 

6 .66 (1.27) 


6.79 

6.67 (1.26) 


Stat Interp 
(Cd = 7) 


7.34 

7.23 ( 1. 25) 


6.87 

6 .75 (1.26) 


6.87 

6.76 (1.25) 


Stat Interp 
(Dmpd Cos) 


7.37 

7.26 (1.28) 


6.91 

6 .79 (1.28) 


6.91 

6 .79 (1 . 28) 


Thin Pi Spl 
(m = 2) 


7.12 

7.00 (1.30) 


6.59 

6 .45 (1.33) 


6.60 

6.46 (1.32) 


Lapl Sm Spl 
(m = 3) 


10.54 

10.40 ( 1. 73) 


10.25 

10 . 10 (1 . 73) 


10.25 

10.10 (1.73) 


Frnke/Gr dn 
(3 Pass) 


12.02 

1 1.72 (2.65) 


11 .75 

11 .45 (2.65) 


11.76 

1 1.45 (2.65) 


PseudoBa rnes * 
(2 Pass) 


9.28 

9. 10 (1.83) 


8.87 

8.68 (1.82) 


8.87 

8.68 (1.8 2) 


PseudoBa rnes ' 
(Lccal) 


8.20 

8.06 (1.51) 


7.70 

7.55 (1. 50) 


7.71 

7.57 (1. 50) 



TABLE 3 



Analysis 
RMSE (StCsv) 



Bsl Bicub 



6.09 

5.98 ( 1. 19) 
6.55 

6.44 ( 1. 18) 

8,88 

8.69 ( 1.82) 
7.96 

7.81 (1.56) 
6.79 

6 .68 ( 1. 26 ) 

6.87 

6.76 ( 1.25) 
6.91 

6.79 ( 1.28) 
6.60 

6.47 ( 1.32) 
10.25 

10.11 (1.73) 
11.76 

11 .45 (2. 65) 

8.88 

8.69 ( 1.82) 
7.72 

7.57(1. jO) 
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E = 500, Theta = 100, Delth = 13.775 Entries: RMSE 
rg = 3C, ro = 10 Mean 
Number of realizations = 100 

13X9 grid of 2.5 degrees, 36 Observation points 



Gr id-to-obs : 


EW Linear 


Bicub Spl 


PW Bicub 


Obs-to-grid 








Opt Inteip 
(Cd = 10) 


6.88 

6.72 ( 1. 49) 


6.19 

6.05 (1 .34) 


6.20 

6.06 (1.34) 


Local 01 
(Cd = 10) 


7.40 

7.24 (1.51 ) 


6.76 

6 .62 (1.37) 


6.77 

6.63 (1.37) 


Bar res * 
2-Eass 


9.73 

9.53 (1.96) 


9.28 

9.09 (1.88) 


9.29 

9.10 (1.88) 


Barnes ' 
(Local) 


8.29 

8.16 ( 1. 48) 


7 .79 

7.67 (1.36) 


7.77 

7.66 (1.36) 


Stat Interp 
(Cd = 14) 


7.71 

7.56 (1.50) 


7 .08 

6.95 (1.35) 


7.09 

6.96 (1.35) 


Stat Interp 
(Cd = 7) 


7.54 

7.39 (1.50) 


6.96 

6.82 (1.36) 


6.96 

6.83 (1.36) 


Stat Interp 
(Diipd Cos) 


7.37 

7.26 (1.28) 


6.91 

6 .79 (1.28) 


6.91 

6.79 (1. 28) 


Thin PI Spl 

(m = 2) 


7.45 

7.29 ( 1. 51) 


6.80 

6.66 (1.37) 


6.81 

6.67 (1.37) 


Lapl Sm Spl 
(m = 3) 








Frnke/Gr dn 
(3 Eass) 








PseudoBa rnes ' 
(2 Eass) 


9.75 

9.55 (1.97) 


9.28 

9.09 (1.88) 




PseudoBa rnes * 
(Local) 









Analysis 

EMSE(StDev) 



Bsl Bicub 



6.21 

6.06 ( 1. 34) 



6.91 

6.79 (1.28) 



TABLE 4 



E = 500, Theta = 100, Celth = 0 Entries: RKSE 
rg - 20, ic = 10 Mean 
Number of realizations = 100 

13X9 grid of 2.5 degrees, 36 Observation points 



Grid-to-cbs : 


EW Linear 


Bicub Spl 


PW Bicub 


Ob£-to-giid 








Opt Inteip 
(Cd = 10) 


6.23 

6. 10 ( 1. 28) 


5.75 

5.62 (1. 22) 


5.76 

5.63 (1.22) 


Local 01 
(Cd = 10) 


6.54 

6.4 1 (1. 28) 


6 . 10 

5.97 (1.22) 


6.10 

5.98 (1.22) 


Barnes ' 
2- Pass 


7.18 

7.03 (1.47) 


6.85 

6 .71 (1.47) 


6.86 

6.71 (1.4 1) 


Barnes ' 
(Lccal) 


7. 19 

7.08 ( 1. 22) 


6 .77 

6.68(1.11) 


6.76 

6.67 (1.10) 


stat Interp 
(Cd = 14) 


6.70 

6.58 ( 1. 27) 


6.30 

6.19(1.19) 


6.31 

6.19 (1.19) 


stat Interp 
(Cd = 7) 


6.71 

6.58 (1. 30) 


6 .24 

6.12 (1.25) 


6.25 

6 . 12 (1. 24) 


stat Interp 
(Dmpd Cos) 


6.78 

6.65 (1.33) 


6.31 

6 . 18 (1.28) 


6.32 

6.18 (1. 28) 


Thin PI Spl 
(m = 2) 


6.66 

6.52 (1. 34) 


6 . 20 

6 .06 (1.29) 


6.20 

6.07 (1.29) 


Lapl Sm Spl 
(m = 3) 


1 1.05 

10.84 (2. 16) 


10.71 

10.50 (2. 12) 


10.71 

10 .50 (2. 12) 


Frnke/Gr dn 
(3 Pass) 


6.72 

8.54 (1.78) 


8.56 

8 .39 (1.70) 


8.56 

8.39 (1. 70) 


Pseudo Barnes ' 
(2 Eass) 


7.19 

7.03 ( 1. 47) 


6.85 

6 . 85 (1.41) 


6.86 

6.71 (1.4 1) 


Pseudo Barnes ' 
(Lccal) 


6.78 

6.66 ( 1. 28) 


6 .37 

6 .25 (1.24) 


6.38 

6.26 (1.23) 



TABLE 5 



Ar clysi 
F.MSE (St 



Bsl Bicub 
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O ^0 



E = 500, Ibeta = 100, Eelth = 0 Entries: RMSE 
ig = 3C, ro = 5 Mean 
Kumber of realizations = 100 

13X9 grid of 2.5 degrees, 36 Observation points 



Grid-to-cfcs ; 


PW Linear 


Bicut Spl 


PW Bicub 


Obs-to-giid 








Opt Inteip 
(Cd = 10) 


4.57 

4.50 (0. 83) 


3 .76 

3.70 (0.69) 


3.77 

3.70 (0.69) 


Local 01 
(Cd = 10) 


5.03 

4.95 (0.91 ) 


4.26 

4 . 19 (0.77) 


4.27 

4.20 (0.77) 


Ba rnes * 
2-Eass 


8.88 

8.69 (1.85) 


8.48 

8 .29 (1.76) 


8.49 

8.30 (1.77) 


Barnes * 
(Lccal) 


6.39 

6.30 ( 1. 10) 


5.90 

5.79(1.11) 


5. 89 

5.78 (1.11) 


St at Inter p 
(Cd = 14) 


5.26 

5. 17 (0.92) 


4.57 

4.49 (0.83) 


4.57 

4 .50 (0.83) 


Stat Interp 
(Cd = 7) 


5.02 

4.95 (0. 82) 


4.28 

4.21 (0.74) 


4.28 

4.21 (0. 74) 


Stat Interp 
(Dmpd Cos) 


4.96 

4.89 (0.83) 


4.22 

4 . 16 (0.75) 


4.22 

4.16 (0. 75) 


Thin Pi Spl 
(m = 2) 


4.92 

4.85 (0.84) 


4.15 

4 .09 (0.72) 


4.16 

4.10 (0.72) 


Lapl Sm Spl 
(m = 3) 


6.08 

5.98 (1.05) 


5.43 

5.33 (1.06) 


5.43 

5.33 (1.06) 


Frnke/Gr dn 
(3 Pass) 


1 1.82 

1 1.55 (2.55) 


11.61 

11.34 (2. 51) 


11.61 

11.34 (2.50) 


Pseudo Barnes’ 
(2 Pass) 


8« 8 9 

8.70 ( 1. 85) 


8 .48 

8.29 (1.76) 


8.49 

8.30 (1.77) 


PseudoBarnes' 

(Local) 


7.30 

7. 15 ( 1. 47) 


6.75 

6 .60 (1. 40) 


6.76 

6.62 (1.41) 



TABLE 6 



Analysis 
BMSE (StDev) 



Bsl Bicub 
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P = 500, Iheta = KANDCf!, Delth = FANDOM, Entries: BMSE 

rg = 30, ro = 10 Mean 

Number cf realizations = 100 

13X9 grid of 2.5 degrees, 36 Observation points 



Gr id-to-cbs : 


PW Linear 


Bicub Spl 


PW Bicub 


Obs-to-grid 








Opt Interp 
(Cd = 1 C) 


7.00 

6.85 (1.44) 




6.40 

6.27 (1.30) 


Local 01 
(Cd = 10) 


7.48 

7.33 ( 1. 49) 




6.90 

6.90 (1.36) 


Barnes ' 
2-Pass 


9.92 

9.70 ( 2. 07) 




9.52 

9.33 (1.89) 


Barnes' 

(Local) 


8.34 

8.21 (1.46) 




7.89 

7.76 (1.39) 


stat Interp 
(Cd = 14) 


7 9 7 

l\ll (1.56) 




7.41 

7.27 (1.43) 


Stat Interp 
(Cd = 7) 


7.56 

7.41 ( 1.51) 




7.01 

6.88 (1.36) 


Stat Interp 
(Dmpd Ccs) 


7.58 

7.43( 1. 53) 




7.04 

6.90 (1.38) 


Thin Pi £pl 
(m = 2) 


7.63 

7.47 ( 1. 58) 




7.06 

6. 92 (1.42) 



Lapl Sm £pl 
(m = 3) 



Fr nke/Grdn 
(3 Pass) 

Pse udoBa rnes * 
(2 Pass) 

PseudoBarnes • 
(Local) 



TABLE 7 



An aly sis 
BMSE (StDev) 



Bsl Bicub 
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E = 500, Theta = BANDCM, Delth = RANDOM, Entries: 
ig = 30, ic = 10, rg (lie) = 20 
Number of realizations = 100 

13X9 grid of 2.5 degrees, 36 Observation points 



RMSE 

Mean 



Gr id-t o-cbs : 


PW Linear 


Obs-to-grid 




Opt Interp 
(Cd = 10) 


7.14 

6.97 (1.52) 


Local 01 
(Cd = 10) 


7.56 

7.39 ( 1.62) 


Barnes ' 
2- Pass 


9.64 

9.42(2. 04) 


Barnes* 

(Local) 


8.30 

7.04 (1.43) 


stat Interp 
(Cd = 14) 




Stat Interp 
(Cd = 7) 




Stat Interp 
(Dnpd Ccs) 




Thin Pi Spl 
(m = 2) 


7.19 

7.04 ( 1. 45) 


Lapl Sm Spl 
(m = 3) 




Frnke/Grdn 
(3 Pass) 




PseudoBaines* 
(2 Pass) 




PseudoBaines • 
(Lccal) 





Bicub Spl PW Bicub 



6.47 

6.32 (1.35) 
6.92 

6.77 (1.41) 

g 21 

9.01 (1. 9C) 
7.73 

7.62 (1.33) 



6.52 

6.39 (1.26) 



TABLE 8 



Analysis 
RMSE (StDev) 



Bsl Bicub 
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E = 500, Ibeta = RANDOM, Delth = RANDOM, Entries: RMSE 
rg = 2C, rc = 10, rg (lie) = 30 Mean 
Number of realizations = 100 

13X9 grid of 2.5 degrees, 36 Observation points 



Gr id-to-cbs : 
Obs-to-grid 


EW Linear 


Bicub Spl 


PW Bicub 


Opt Interp 


6.32 




5.72 


(Cd = 10) 


6.17 (1.36) 




5.60 (1 .2 1) 


Local 01 


6.60 




6.05 


(Cd = 10) 


6.47 ( 1. 33) 




5. 94 (1. 16) 


Barnes ' 


7.10 




6.69 


2- Eass 


6.95 ( 1. 45) 




6.55 (1.39) 


Barnes* 


7.09 




6.54 


(Local) 

Scat Interp 
(Cd = 14) 

stat Interp 
(Cd = 7) 

stat Interp 
(Dipd Ccs) 


6.98 (1.22) 




6.45 (1. 08) 


Thin PI £pl 


6.31 




5.75 


(m = 2) 

Lapl Sm £pl 
(m = 3) 

Fr nke/Gr dn 
(3 Eass) 

PseudoBarnes ' 
(2 Eass) 

PseudoBarnes ' 
(Local) 


6.17( 1.30) 




5.63 (1.16) 



TABLE 9 



Analysis 
RMSE (StDev) 



Bsl Bicub 
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E = 500, Theta = RANDCM, Delth = RANDOM, Entries: RMSE 
rg = 3C, rc = 10, ro (lie) = 5 Mean 
Number of realizations = 100 

13X9 grid of 2.5 degrees, 36 Observation points 



Grid-to-cbs: 


PW Linear 


Obs-to-grid 




Opt Interp 
(Cd = 10) 


7.51 

7.37 (1.43) 


Local 01 
(Cd = 10) 


7.96 

7.85( 1. 36) 


Barnes ' 
2- Pass 


9.81 

9.60(2.01) 


Barnes' 

(Local) 


8.47 

8.36 (1.32) 


Stat Interp 
(Cd = 14) 




Stat Interp 
(Cd = 7) 




Stat Interp 
(Dmpd CCS) 




Thin PI Spl 
(m = 2) 


7.65 

7.51 ( 1. 43) 


Lapl Sm Spl 
(m = 3) 




Frn)?e/Gr dn 
(3 Pass) 




PseudoBarnes ' 
(2 Pass) 




PseudoBarnes ' 
(Local) 





Bicub Spl PW Bicub 



6.88 

6.74 (1.38) 
7.37 

7.26 (1.27) 
9.28 

9.06 (1.99) 
7.91 

7.80 (1.31) 



6.94 

6.80 (1.39) 



TABLE 10 



Analysis 
RMSE (StDev) 



Bsl Bicub 
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E = 500, 


Theta = RANDOM, Delth = 


RANDOM, 


Entries: RMSE 


ig = 30, 


rc = 5 , ro (lie) =10 




Mean 


Number of 


realizations = 100 






13X9 grid 


of 2.5 degrees, 36 Observation 


points 


Gr id-to-cbs 


: FW Linear Bicub 


Spl 


PW Bicub 


Obs-to-grid 








Opt interp 


5. 18 




4.19 


(Cd = 1C) 


5.10 (0.91) 




4.13 (0.7 1) 


Local 01 


5.66 




4. 76 


(Cd = 10) 


5.56 ( 1. 04) 




4.6 9 (0.81) 


Barnes * 


9.09 




8.62 


2-Eass 


8.93( 1. 69) 




8.47 (1.61) 


Baxres' 


6.71 




6.06 


(Local) 


6.63 (1.07) 




5.97 (1.04) 


Stat Interp 
(Cd = 14) 








Stat Interp 
(Cd = 7) 








stat Interp 
(Diipd Ccs) 








Thin Pi £pl 


6.13 




5.34 


(m = 2) 


6.02( 1. 15) 




5.25 (0.95) 


Lapl Sm £pl 
(m = 3) 








Frnke/Grdn 
(3 lass) 








PseudoBaine 
( 2 Pas s ) 


g t 






PseudoBaine, 

(Local) 


c 1 







TABLE 11 



Analysis 
RMSE (StDev) 



Bsl Bicah 
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P = 50C, Theta = EANDCK, Delth = BANDOM, Entries: 

rg = 5, rc = 10 

Number cf realizations = 100 

13X9 grid of 2.5 degrees, 36 Observation points 



BMSE 

Mean 



Grid-to-cbs : 


PW Linear 


Bicub Spl 


PW Bicub 


Obs-to-grid 








Opt Interp 
(Cd = 1 C) 


3.58 

3.44 (1.01) 




3.32 

3.18 (0. 94) 


Local 01 
(C d = 10) 


3 77 

3*63 ( 1. 02) 




3.53 

3.40 (0. 9769 


Barnes • 
2-Eass 


4.44 

4.32( 1. 0 1) 




3 98 

3.87 (0. 90) 


Barnes' 

(Local) 


6.26 

6.15 (1.16) 




5.73 

5.64 (1.00) 


Stat Interp 
(Cd = 14) 


3.61 

3.46 (1.03) 




3.41 

3.27 (0.96) 


Stat Interp 
(Cd = 7) 


3.65 

3.50(1.03) 




3.39 

3.25 (0.97) 


Stat Interp 
(Dnpd Ccs) 


3.78 

3.62( 1.07) 




3.51 

3.36 (1.02 


Thin PI £pl 
(m = 2) 


4.00 

3.84 ( 1. 13) 




3.85 

3.70 (1.07) 


Lapl Sm £pl 
(m = 3) 








Frnke/Gr dn 
(3 Pass) 








PseudoBa rnes * 
(2 Pass) 








PseudoBaines ' 
(Local) 




TABLE 12 





Analysis 
BMSE (StDev) 



Bsl Bicub 
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E = 500, Iheta = RANDOM, Delth = RANDOM, Entries: RMSE 

ig = 30, rc = 0 Mean 

Number of realizations = 100 

13X9 grid of 2.5 degrees, 36 Observation points 



Grid-to-cbs: EW Linear Bicut Spl PW Bicub 



Obs-to- grid 



Opt Interp 


6.06 


0.74 


(Cd = 10) 


5.84 ( 1. 62) 


0.69 (0.27) 


Local 01 


3.86 


2.40 


(Cd = 1 0) 


3.79 (0.76 ) 


2.2. 30(0.69 


Barnes ' 


9.38 


8.96 


2-Pass 


9.16 (2.05) 


8.74 (1.98) 


Barnes ' 


6.27 


5.63 


(Lccal) 


6. 15( 1. 21) 


5.49 (1.23) 


Stat Interp 


7.62 


1.23 


(Cd = 14) 


7.10 (2.74) 


1 . 16 (0. 40) 


Stat Interp 


3.70 


1.29 


(Cd = 7) 


3.64 (0.68) 


1 .20 (0. 48) 


Stat Interp 


4.48 


1 .03 


(Dmpd Cos) 


4.37 (0. 96 ) 


0.97 (0.36 


Thir. PI Spl 


3.74 


2.18 


(m = 2) 


3.66 (0.78) 


2.02 (0.83) 



Lafl Sm Spl 
(m = 3) 

Err ke/Gr dn 
(3 Pass) 

Pseudo Earnss ' 
(2 Pass) 

Pseudo Earnes’ 
(Lccal) 



TABLE 13 



An aly si 
RMSE (St 



Bsl Bicub 
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(/) o 



Entries: 



RMSE Analysis 
Mean RMSE(StDev) 



E = 500, Theta = 100, Celth = 0 
ig = 30, lo = 10 
Number of realizations = 100 
11 BY 8 grid of 5 degrees, 67 Observation points 



Gr id-to-obs ; 


PW Linear 


Bicub Spl 


PW Bicub 


Bsl Bicub 


Obs-to-giid 










Opt Interp 
(Cd = 10) 


12.84 

12.74 (1.62) 


7.62 

7.53 (1.19) 


7.93 

7.84 (1.20) 


8.11 

8.02(1. 19) 


Local 01 
(Cd = 10) 


13.33 

13.22 (1.73) 


8.44 

8 .33 (1. 33) 


8.74 

8.63 (1.35) 


8.92 

8.82 ( 1. 32) 


Barnes* 

2-Pass 


14.33 

14.21 (1.85) 


10 .62 

10.49 (1.70) 


10.82 

10.68 (1.70) 


10.95 

10.81 (1.71) 


Barnes * 
(Local) 


14.00 

13.91 ( 1. 55) 


8.82 

8.71 (1.36) 


<5 12 

9.79 (1.35) 


9.40 

9.31 (1. 33) 


Stat Interp 
(Cd = 14) 


13.75 

13.25( 1. 57) 


8.80 

8.70 (1 .33) 


9.02 

8.92 (1.31) 


9.20 

9.10 (1. 32) 


Stat Interp 
(Cd = 7) 


13.44 

13.35 (1. 63) 


8.31 

8.23(1.18) 


8.62 

8.53 (1.21) 


8.79 

8.70( 1.19) 


Stat Interp 
(Dmpd Cos) 


13.57 

13.47 (1.65) 


8.47 

8 .38 (1. 25) 


8.78 

8.68 (1.27) 


8.95 

8.86 (1.25) 


Thin PI Spl 
(m = 2) 


13.17 

13.07 (1.59) 


8.08 

7.99 (1 . 17) 


8.36 

8.27 (1.19) 


8.47 

8.39 (1. 19) 


Lapl Sm Spl 
(m = 3) 


17.12 

17.01 ( 1. 86) 


11 .87 

11 .76 (1.60) 


12.08 

1 1 . 97 (1.57) 


12.16 

12.05 (1.58) 


Frnke/Gr dn 
(3 Eass) 


17.29 

17. 14 (2. 29) 


15.23 

15.04 (2.35) 


15.31 

15.13 (2.34) 


15.36 

15.17 (2.37) 


PseudoBarnes ' 
(2 Eass) 


14.27 

14.14 (1.87) 


10 .62 

10.49 (1.70) 


10.82 

10.69 (1.70) 


10.95 

10.82 (1.71) 


PseudoBarnes * 
(Lccal) 


13.89 

13.77 (1. 82) 


9.73 

9 .60 (1. 58) 


10.00 

9.87 (1.58) 


10.17 

10.05 (1.57) 



TABLE 14 
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Estimated end (empirical) RMS errors for statistical methods 



ig: 
ro : 

method 


5 

10 


20 

10 


30 
1 0 


30 

5 


30 

0 


Opt Interp 


3.30 


5.62 


6.29 


3. 78 


0.65 


(cd = 10) 


(3. 32) 


(5.76) 


(6. 40) 


(3.77) 


(0. 74) 


Local 01 


3. 55 


5.96 


6. 80 


4.39 


2. 21 


(cd = 10) 


(3.5 3) 


(6.10) 


(6.90) 


(4. 27) 


(2.40) 


Stat Inteip 


2.90 


4.60 


5. 07 


2. 95 


0.11 


(cd = 14) 


(3.41) 


(6.31) 


(7.41) 


(4.57) 


(1. 23) 


Stat Inteip 


3. 76 


7. 14 


8. 26 


5.46 


2. 90 


(cd = 7) 


(3.39) 


(6.25) 


(7.01) 


(4. 28) 


(1. 29) 


Stat Interp 


3.82 


7.00 


7. 94 


5. 03 


2. C6 


(Dmpd Cos) 


(3.51) 


(6.32) 


(7.04) 


(4.22) 


(1.03) 



TABLE 15 



2.5 DEGREE GRID fiND OBSERVATION LOCATIONS 




LO 

CD 



LO 

Cv. 




LD 

OO 



LO 

CD 



LD 

O 



o 

LO 



N „09 



N „0S 



N 



o 



ID 

CM 




' ID 

N gOC r 52 

N „0S 
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Figure 



5 DEGREE GRID AND OBSERVATION LOCATIONS 




Ln 

CD 




LO 



LO 

00 



in 
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N 



. 01 ^ 
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105 



500 MB CONTOURS, THETR = 100, DELTH 



CD 

11 







CO 
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Figure 



SKILL (= 1 - RMSfl/RMS OBS) 

- 0.50 - 0.25 0.00 0.25 0.50 



T3: RG = 30, RO = 10, 13X9 2.5 OLGREE GRID 



LEGEND 

01 

LOCRL 01 
BRRNES 

LOORL BRRNES 

SI-CD=7 

TPS 



1 \ 1 I I 

0.0 50.0 100.0 150.0 200.0 250.0 

COST (K-OPS) 




Figure 4 
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SKILL (= 1 - RMSfl/RMS OBS) 

- 0.50 - 0.25 0.00 0.25 0.50 



T5: RG = 20, RO = 10, 13X9 2.5 DEGREE GRID 



LEGEND 
□ = 01 

o = LOGRL 01 
A = BRRNES 
+ = LOCRL BRRNES 
X = SI-CD=7 
o = TPS 



I I I I I 

0.0 50.0 100.0 150.0 200.0 250.0 

COST (K-OPS) 



Q 



P 
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T 

4> 



I 



Figure 5 
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SKILL (= 1 - RMSR/RMS OBS) 

- 0.50 - 0.25 0.00 0.25 0.50 



T6: RG = 30, RO = 5, 13X9 2.5 DEGREL GRID 



□ 

Q X 
[] 

m X 



' ’ " I 

0.0 50.0 



LEGEND 
□ = 01 

o = LOCRL 01 
A = BRRNES 
+ = LOCRL BRRNES 
X = SI-CD=7 
o = TPS 



— ^ 1 1 j 

100.0 150.0 200.0 250.0 

COST (K-OPS) 

Figure 6 
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SKILL (- 1 - RMSR/RMS OBS) 

-0.50 -0.25 0.00 0.25 0.50 



T7: RG = 30, RO = 10, 13X9 2.5 DEGREE GRID 



CP 



Q 



6 



<j> 

4 



+ 



+ 



LEGEND 
n = 01 

o = LOGRL 01 
A = BRRNES 
+ = LOCRL BRRNES 
X = SI-CD=7 
o = TPS 



0.0 



50.0 



100.0 

COST 



150.0 

(K-OPS) 



200.0 



250.0 



Figure 7 
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SKILL (= 1 - RMSn/RMS OBS) 

- 0.50 - 0.25 0.00 0.25 0.50 



TH: RG = 30, RO = 10, 11X8 5 DEGREE GRID 



Q 



6 



LEGEND 
□ = 01 

o = LOCAL 01 
A = BARNES 
+ = LOCAL BARNES 
X = SI-CD=7 
o = TPS 



A 



-n — 

200.0 



0.0 



50.0 100.0 150.0 

COST (K-OPS) 



250.0 



Figure 8 
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RMS flNRLYSIS ERROR RMS ANALYSIS ERROR 

0.0 5.0 10.0 0.0 5.0 10.0 



RO = 10, 13X9 2.5 DEGREE GRID 




Figure 9 



RG = 3D, 13X9 2.5 DEGREE GRID 




Figure 10 
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RMS ANALYSIS ERROR 

5.0 10.0 



I'l 




Figure 11 
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RMS RNRLYSIS ERROR RMS RNRLYSIS ERROR 



PERPORMRNCE DEGRRDRTION, RG = 30, RO = 5 



o 




figure 12 



PEREORMRNCE DEGRRDRTION, RG = 20, RO = 10 





==4 
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LEGEND 




□ = 01 




o = LOCRL 01 




A = BRRNES 




+ = LOCRL BRRNES 




o = TPS 



NOLIE RGLIE=30 



Figure 13 
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